Evidence for Superfluidity of Ultracold Fermions in an Optical Lattice 
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The study of superfluid fermion pairs in a periodic potential has important ramifications for un- 
derstanding superconductivity in crystalline materials. Using cold atomic gases, various condensed 
matter models can be studied in a highly controllable environment. Weakly repulsive fermions in an 
optical lattice could undergo d-wave pairing [l| at low temperatures, a possible mechanism for high 
temperature superconductivity in the cuprates |2(. The lattice potential could also strongly increase 
the critical temperature for s-wave superfluidity. Recent experimental advances in the bulk include 
the observation of fermion pair condensates and high-temperature superfluidity 0, El IE El El 13 • 
Experiments with fermions [E [IE El and bosonic bound pairs [IE in optical lattices have been 
reported, but have not yet addressed superfluid behavior. Here we show that when a condensate 
of fermionic atom pairs was released from an optical lattice, distinct interference peaks appear, 
implying long range order, a property of a superfluid. Conceptually, this implies that strong s-wave 
pairing and superfluidity have now been established in a lattice potential, where the transport of 
atoms occurs by quantum mechanical tunneling and not by simple propagation. These observations 
were made for unitarity limited interactions on both sides of a Feshbach resonance. For larger lat- 
tice depths, the coherence was lost in a reversible manner, possibly due to a superfluid to insulator 
transition. Such strongly interacting fermions in an optical lattice can be used to study a new class 
of Hamiltonians with interband and atom-molecule couplings [lj . 



Previous experiments demonstrating long-range phase 
coherence in Bose-Einstein condensates (BECs) and in 
fermion superfluids used ballistic expansion to observe 
interference of two independent condensates [l5|, vortex 
lattices [E EH 113 or interference peaks after release from 
an optical lattice However, for strongly inter- 

acting fermions, elastic collisions can change the momen- 
tum distribution and wash out interference peaks. For 
an initially superfluid cloud, such dissipative dynamics 
corresponds to superfluid flow faster than the critical ve- 
locity. Consistent with this expectation, expansion of 
the strongly interacting Fermi gas from an optical lattice 
yielded a diffuse cloud which exhibited no signs of inter- 
ference (Fig. 1). This issue was addressed by using a 
magnetic field ramp which quickly increased the detun- 
ing from a Feshbach resonance, taking the system out 
of the strongly interacting regime and enforcing ballis- 
tic expansion. In previous studies of strongly interacting 
Fermi gases magnetic field sweeps were applied to pre- 
vent fermionpairs above the Feshbach resonance from 
dissociating [E, 0, H(j- I n contrast, our experiment re- 
quired a magnetic field sweep both above and below the 
Feshbach resonance to avoid elastic collisions. 

Our experiments used a balanced mixture of 
6 Lifermions in the two lowest hyperfine states. Evapo- 
rative cooling produced a nearly pure fermion pair con- 
densate which was adiabatically loaded into a three- 
dimensional optical lattice. A broad Feshbach resonance 
centered at 834 G enabled tuning of the interatomic in- 
teractions over a wide range. On resonance, a bound 
molecular state becomes degenerate with the open atomic 
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scattering channel, leading to a divergence in the scatter- 
ing length a. Here we explore the region of strong inter- 
actions, also known as the BEC-BCS (Bardeen-Cooper- 
Schrieffer) crossover, where the magnitude of the inter- 
action parameter \kpa\ is greater than unity, and kp is 
defined as the peak Fermi wavevector of a two-component 
non-interacting mixture of 6 Li atoms. Throughout this 
regime, pairing occurs as a result of many-body interac- 
tions. Below resonance, for strong interactions, the bare 
two-body state has a bond length larger than the inter- 
atomic spacing and is irrelevant. In a lattice, atom pairs 
above and below the resonance can be confined to one 
lattice site [ll| , and crossover physics may require an oc- 
cupation larger than or equal to one. 

The peak pair filling factor of the lattice was about 
unity. At this density in the bulk, the fermion pair size 
is on the order of 1/kp = 170 nm, comparable to the lat- 
tice spacing of 532 nm. To probe the momentum distri- 
bution, we ramped the magnetic field out of the strongly 
interacting regime as fast as technically possible (~ 150 
/is) and then turned off the confining potential. Absorp- 
tion images taken after 6.5 ms of expansion reveal sharp 
peaks at the reciprocal lattice vectors — the signature of 
long-range coherence and superfluidity. 

We observed such interference peaks at magnetic fields 
both above and below the Feshbach resonance (Fig. 2). 
The six first order diffracted peaks are clearly visible 
around the zero momentum fraction and their positions 
correspond to the expected momentum quanta of 2hkL 
carried by molecules of mass 2m, where ki, is the lat- 
tice wavevector. At high magnetic fields (Fig. 2d) the 
visibility of the interference peaks decreased and some 
additional heating was observed. This degradation could 
be due to a higher fraction of thermal atoms as we ap- 
proached the BCS limit, but was not studied in detail. 
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FIG. 1: The schematic shows the time sequence of the 
magnetic field ramp used throughout this Letter. A one- 
dimensional optical standing wave was pulsed onto the super- 
fluid at different values B p (indicated by arrows (a) 822 G, (b) 
749 G and (c) 665 G) of the magnetic field during expansion, 
creating particles at twice the photon recoil |2g|. Absorption 
images taken at the time marked with the cross show distinct 
momentum peaks only at magnetic fields B p < 750G, corre- 
sponding to kFd < 1 . At higher magnetic fields, the peaks 
blurred into a broad diffuse cloud as a result of the larger 
collision cross-section. 






FIG. 2: (a), The orientation of the reciprocal lattice, also 
with respect to the imaging light, b-d, Interference peaks 
are observed for magnetic fields of (b) 822 G, (c) 867 G and 
(d) 917 G. The lattice depth for all images is 5 E r and each 
image is the average of 3 shots. The field of view is 1 mm 
by 1 mm. Density profiles through the vertical interference 
peaks are shown for each image. 



The narrow interference peaks clearly reveal the pres- 
ence of a macroscopic wavefunction possessing long range 
phase coherence. The separation between the interfer- 
ence peaks relative to their width gives an estimate of 
the coherence length of ~ 10 lattice sites. This estimate 
is a lower bound, because effects of finite resolution and 
mechanisms of residual broadening have been neglected. 
With unity occupation, and in the absence of any dis- 
cernible background at magnetic fields near the Feshbach 
resonance, this implies a minimum phase space density of 
10 3 , and shows that our samples are deep in the quantum- 
degenerate regime. In previous studies of ultracold Bose 
and Fermi gases, the appearance of a condensate frac- 
tion and long range phase coherence was shown to oc- 
cur concurrently with the possibility to excite superfluid 
flow H, EH 0, E[ . Superfluid hydrodynamics is usually 
regarded as the direct proof for superfluidity. However, 
all reports of superfluidity of bosons in three-dimensional 
optical lattices have relied solely on observations of sharp 
interference peaks and inferred superfluidity from the es- 
tablished connection between long-range coherence and 
superfluidity 0,[2^. Similarly, our observations directly 
show long-range coherence and indirectly show superflu- 
idity of fermion pairs in an optical lattice. 

For deep lattices, breakdown of superfluid behavior 
has been observed pj| 12*1 fc> r weakly interacting Bose- 
Einstein condensates. This phase transition to the Mott- 



insulator state occurs when on-site interactions start 
to suppress atom number fluctuations and the system 
transitions from a delocalized superfluid described by 
a macroscopic wavefunction, to a product of Wannier 
states tightly localized at each lattice site. Experimen- 
tally, this is manifested as a smearing of the distinct 2hkL 
interference peaks. 

Fig. 3 shows the evolution of a strongly interact- 
ing fermionic superfluid when the lattice depth was 
increased. The interference peaks became more pro- 
nounced initially, due to increased modulation of the 
wavefunction. The interference peaks began to smear 
out, rapidly giving way to a featureless cloud, beyond a 
critical lattice depth V c « 6 E r where E r = h 2 k\/Am = 
hx 15 kHz is the recoil energy. This indicates that all 
phase coherence had been lost. Upon subsequent ramp 
down of the lattice, interference peaks became visible 
again (Fig. 3h), demonstrating reversibility of the lat- 
tice ramp. 

We repeated this sequence for a wide range of initial 
magnetic fields both above and below the resonance and 
observed the same marked change in the interference pat- 
tern. Fig. 4 displays the peak optical density of the 
interference peaks for different lattice depths at repre- 
sentative fields. Across all fields, the sharp decrease in 
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FIG. 3: Values of V 's are (a) E r , (b) 2.5 E r , (c) 4 E r 
, (d) 5 E r , (e) 6 E r , (f) 7 £ r , (g) 9 E r , (h) 2.5 £ r . 
(a-g) were taken after an adiabatic ramp up to the final Vb, 
while (h) was taken after first ramping up to 10 E r , before 
ramping down to 2.5 E r . 
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FIG. 4: Values of magnetic fields are 842 G (filled circles), 
892 G (open squares) and 942 G (filled triangles). Peak op- 
tical densities were estimated from fits to the peaks, includ- 
ing background subtraction. The inset (color online) shows 
a sample density profile of the central and one pair of inter- 
ference peaks (black dotted line), with a bimodal fit to one 
side peak (red solid line). Each point is the average of three 
different images with six interference peaks per image. Error 
bars show s. d. 



peak optical density occurred between 5 and 6 E r . Fur- 
ther increase of the magnetic field resulted in decreasing 
overall visibility, until interference peaks could no longer 
be observed regardless of lattice depth. 

The loss of phase coherence with increasing lattice 
depth is consistent with the qualitative description of 
the superfluid to Mott-insulator transition. However, 
the usual single-band description is no longer applicable, 
since in the strong-coupling regime the on-site interaction 
strength should be comparable to the band gap Hlj, where 
uo is the onsite trap frequency. Furthermore, Pauli block- 
ing forbids the multiple occupation of the lowest state 
of an individual lattice site by identical fermions, and 



160- 



140- 



03 
QJ 
Cl 



c 

QJ 

u 



120- 



100- 



80 - 



60 



8 8 



o o - 
o o o o 
Q o 



X 



t ♦ 
♦ 



♦ ♦ 



200 



400 



600 



—I 

800 1 000 



Ramp down time (|is) 



FIG. 5: The width of the central peak is used as a measure of 
phase coherence after an adiabatic ramp up to 8 E r followed 
by a fast ramp down to 2.5 E r at a fixed magnetic field of 822 
G. Filled circles were extracted with the use of a gaussian fit, 
and diamonds with a bimodal fit. Also plotted for compari- 
son is the gaussian width of the central peak for a dephased 
sample, in which a field gradient was applied during the ramp 
up of the lattice (open circles). All points were taken for 6.5 
ms time-of- flight. 



modification of the single particle tunneling rate is ex- 
pected due to virtual pair breaking transitions 14]. One 
may still be tempted to use the standard bosonic Hub- 
bard model and estimate the critical lattice depth V c for 
an assumed value of onsite interaction energy U = Huj 
and non-interacting, single particle tunneling J, but the 
obtained V c ~ 3 E r is significantly smaller than our ob- 
servation, which is in turn much smaller than the Vr > 10 
E r observed for weakly interacting atomic BECs [l^J, 123 • 
Along with the observed insensitivity of V c to magnetic 
field, this demonstrates that models based on weak in- 
teractions are inadequate. 

Fig. 3h demonstrates the reversibility of the transi- 
tion from a long range coherent state to a state without 
strong coherence. We now study the time scale for this 
recoherence, in analogy with similar measurements per- 
formed across the superfluid to Mott-insulator transition 
in atomic BECs 19]. Fig. 5 shows that phase coher- 
ence was restored on a sub-millisecond time scale, on the 
order of the the single-particle tunneling time of about 
500 /is (for a shallow lattice of 2.5 E r ). When the same 
lattice ramp sequence was applied to a superfluid which 
had been dephased by a magnetic field gradient pjj the 
system did not regain phase coherence on the time scales 
that we probed. Therefore, evaporative cooling is neg- 
ligible during this time. The short recoherence time of 
the condensate is evidence that the system stayed in its 
ground state or at least in a low entropy state when the 
lattice was ramped up. 

Fig. 5 also provides evidence that the system could 
not recohere during the 150 /is magnetic field ramp. In 
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Fig. 3h, the central peak is well fitted by a bimodal dis- 
tribution with a width of 35 /im, in clear contrast to the 
gaussian width of 105 /im obtained from Fig. 5 after 150 
/is. Therefore, we conclude that the observed interfer- 
ence patterns in Fig. 1 reflect the coherence of the cloud 
at the initial magnetic field, in the strongly interacting 
regime. 

We have shown long range phase coherence of fermion 
pairs in an optical lattice in the BEC-BCS crossover re- 
gion by observing sharp interference peaks during bal- 
listic expansion. This indicates that we have achieved 
strong s-wave pairing and superfluidity in a lattice po- 
tential. Further studies will reveal how the pair wave- 
function is affected by confinement [24[ , and whether the 
lattice shifts the BEC-BCS crossover away from the Fesh- 
bach resonance |2a| . The loss of coherence during the lat- 
tice ramp up and the rapid recoherence are characteristic 
of a Mott-insulator. However, definitive proof requires 
a better understanding of the unitarity-limited interac- 
tions in such a Fermi system. Recent theoretical work 
[l4l predicts that strongly interacting fermions in 
an optical lattice feature multi-band couplings and next- 
neighbor interactions and can realize the important t — J 
and magnetic XX Z models of condensed matter theory. 
This demonstrates that such atomic systems are an ideal 
laboratory for the exploration of novel condensed-matter 
physics. 

Methods 

Clouds of superfluid fermion pairs were created in a 
new experimental setup j23, using techniques similar 
to those described elsewhere |8|. In brief, a combination 
of laser cooling and sympathetic cooling of spin polar- 
ized fermions by bosonic 23 Na was followed by a spin 
transfer to create a two component Fermi gas, allowing 
further cooling via direct evaporation of the fermions. 
As the fermions cooled, they formed pairs which Bose- 
condensed. 

Estimates of the scattering length and hence the inter- 
action parameter from magnetic field were obtained using 
a(B) = -1405a (l + 300/(£-834))(l + 0.0004(£-834)) 
|29j |. where B is measured in Gauss, and ao is the Bohr 
radius. The calibration of the magnetic field in our sys- 
tem had an uncertainty of about 5 G. 

Evaporation was performed at a magnetic field of 822 
G, where strong interactions allowed for efficient evapo- 



ration. An estimated average final number of N c± 2 x 10 5 
6 Li pairs and harmonic trapping frequencies of v x ,y,z = 
(270, 340, 200) Hz gave a trap depth of 1.7 /iK and Fermi 
energy E F = ks x 1.4 /iK, where E F = fWJ^N) 1 ^ and W 
is the average trapping frequency. After evaporation, the 
magnetic field was brought to a desired value Bo in 20 ms 
and the condensate allowed to equilibrate for a further 
200 ms. Before ramping to values of Bo on the BCS side, 
we also recompressed the optical trap to (340,440,270) 
Hz and 2.2 /iK depth in 100 ms to accommodate the 
larger Fermi clouds above the resonance Q. 

A three-dimensional optical lattice was formed from 
three optical standing waves, oriented such that the re- 
sulting unit cell has a sheared cubic structure, with one 
axis tilted ~ 20 degrees from the normal for reasons of op- 
tical access (see Fig. la) 23]. The incident laser beams 
were focused down to the condensate with waists of ~ 
90 /im, then retro-reflected and overlapped at the con- 
densate to generate the standing wave potentials. All 
lattice light was derived from a 1064 nm single frequency 
fiber laser, and each beam was detuned tens of MHz with 
respect to the others to eliminate interference between 
different beams. 

The lattice potential was imposed onto the conden- 
sate by adiabatically increasing the intensity of the laser 
beams to a variable final value Vq. The calibration of 
Vo had an uncertainty of about 20 percent. A simple 
linear ramp with a constant rate dVo/dt of 0.5 E r per 
ms was used unless otherwise specified. This satisfies the 
interband adiabaticity condition of dVo/dt « 16E^/h. 

Ballistic expansion for the detection of the different 
momentum components was provided by a magnetic field 
sequence (shown in Fig. 1) that quickly brought the sys- 
tem out of the strongly interacting regime when all con- 
finement was switched off. During the magnetic field 
ramp of about 150 /is, the lattice potential was kept 
on. The first 2 ms of expansion took place at 470 G, 
where the molecules are tightly bound, before the field 
was ramped back up to 730 G in the next 4.5 ms, at 
which the weakly bound molecules strongly absorb light 
near the atomic resonance line and could be observed by 
absorption imaging. The specific magnetic field sequence 
was chosen to minimize collisions within technical capa- 
bilities. 
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